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7^, ' Saturation of small-x gluons in a nucleus, which has the form of transverse 

momentum broadening of projectile gluons in pA collisions in the nuclear rest 
frame, leads to a modification of the parton distribution functions in the beam 



C^ \ compared with pp collisions. The DGLAP driven gluon distribution turns 



out to be suppressed at large x, but significantly enhanced at a; <^ 1. This is 
a high twist effect. In the case of nucleus-nucleus collisions all participating 
[■~~. ■ nucleons on both sides get enriched in gluon density at small x, which leads to 

O \ a further boosting of the saturation scale. We derive reciprocity equations for 

the saturation scales corresponding to a collision of two nuclei. The solution 
of these equations for central collisions of two heavy nuclei demonstrate a 
significant, up to several times, enhancement of QIa-, in AA compared with 
pA collisions. 
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1. Introduction 

The transverse momentum distribution of gluons in nuclei is known to 
be modified compared with a free nucleon. The mean transverse momentum 
squared increases up to a value called saturated scale, Ql^, which depends 
on the nuclear profile. This phenomenon, called color glass condensate [l|, 
is related to parton saturation at small x [2|, and can be also understood in 
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terms of the Landau- Poineranchuk principle [3j as a consequence of coherent 
gluon radiation from multiple interactions in the nucleus [^ . The value of the 
saturation momentum was calculated and compared with data on broadening 
in [^ and has been modeled recently in p, |6|, [tI, |8|, M, llO, lul, 112 ■ 



The saturation scale can be measured as pr-broadening of a parton prop- 
agating through the nucleus in its rest frame [4] , 

Q'sA = ^pI (1) 

Although in leading order both sides of this relation rise linearly with nuclear 



profile Ta [13|, [Ij, [151, [161, Il7[, this dependence slows down by gluon shad- 
owing. These phenomena, broadening and suppression of gluons, are closely 
related, since both result from coherence of gluon radiation in multiple inter- 
actions. Solving the corresponding equation derived in y, one arrives at a 
saturation scale considerably reduced compared to the leading order. The Ta 
dependence is slower than linear, and at very large (unrealistic) nuclear thick- 
nesses the saturation scale saturates, becoming independent of Ta- Notice 
that the solution found in [^ is similar to the result of numerical solution 22 
of the Balitsky-Kovchegov equation 18|, ll9| . Broadening of gluons radiated 



in heavy ion collisions was studied with numerical simulations in j20l . |21 



2. Modification of the beam PDF by a nuclear target 

Due to broadening a nuclear target probes the parton distribution in the 
beam hadron with a higher resolution. Therefore, the effective scale Q^ for 
the beam PDF drifts to a higher value Q^ + Q'^a- -^^ ^^^^ glance this seems 
to contradict casuality, indeed, how can the primordial parton distribution 
in the hadron depend on the interaction which happens later? However, 
there is nothing wrong. The interaction performs a special selection of Fock 
states in the incoming hadron. The same phenomenon happens when one is 
measuring the proton parton distribution in DIS. The proton PDF "knows" 
in advance about the virtuality of the photon which it is going to interact 
with. 

The shift in the scale also can be interpreted as a manifestation of the 
Landau- Pomeranchuk principle [3]: at long coherence times gluon radiation 
(which causes the DGLAP evolution) does not depend on the details of mul- 
tiple interactions, but correlates only with the total momentum transfer, 
q + ApT, which after squaring and averaging over angles results in Q'^ + Ap'^. 



As far as the PDF of the projectile proton has a harder scale in pA 
collisions than in pp, the ratio of parton distributions should fall below one 
at forward and rise above one at backward rapidities. This may look like a 
breakdown of fc^-factorization, however, it is a higher twist effect 

Examples of pA to pp ratios Ra{x, Q^ 
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calculated with MSTW2008 
are shown in Fig. [T] for d-quark and gluon distributions in a hard reaction 
[high-pT, heavy flavor production, etc.). We see that the shift in the hard 





Figure 1: Ratio of parton distribution functions in a reaction characterized by a hard scale 
Q^ = 2, 3, 5, lOGeV^ on a nuclear (A=200) and proton targets. Left panel: ratio of the 
d-quark distributions for the quark saturation momentum Q^^ = 1.2 GeV . Right panel: 
Ratio for gluons with the gluon saturation momentum Ql^ = 2 GeV . 

scale caused by saturation in the nucleus leads to a sizable suppression in 
the projectile parton distribution at large x — > 1 and enhancement at small 
a; ^ 1. We also observe that the magnitude of nuclear modification quickly 
decreases with Q^ confirming that this is a high twist effect. 

Important for what follows is the observation of a considerably increased 
population of small-x partons in the projectile proton in pA compared with 
pp collisions. 



3. Nucleus-nucleus collisions: reciprocity relations 

Notice that in pA collisions the modification of the PDFs of the beam 
and target are not symmetric. Namely, the scale of the PDF of the beam 



proton gets a shift, Q^ =^ Ql^^ = Q"^ + Ql^, while the PDFs of the bound 
nucleons, which do not undergo multiple interactions, remain the same as in 
pp collisions. 

The situation changes in the case of a nucleus-nucleus collision: the bound 
nucleons in both nuclei participate in multiple interactions, therefore the 
scales of PDFs of all of them are modified. However, this modification goes 
beyond the simple shift Q^ ^ Q^ + Q\j^. Indeed, in an AB nuclear collision 
not only the two nucleons (one from A and one from S) participating in the 
hard reaction undergo multiple interactions, but also many other nucleons, 
the so called participants, experience multiple soft interactions. For this 
reason their parton distributions are boosted from the soft scale ^ up to 
the saturation scale fi"^ ^ fi'^ + Qly^rs)' which is usually much larger. Thus, 
the participant nucleons on both sides are boosted to a higher scale and get 
softer PDFs, with larger parton multiplicities at small x. This is illustrated 
on the cartoon in Fig. [2l 





Figure 2: Left: pA collision in which the colliding proton is excited by multiple interactions 
up to a saturated scale Q^^ , what leads to an increased multiplicity of soft gluons in the 
incoming proton. Right: nuclear collision in which participating nucleons on both sides 
are boosted to the saturation scales, Q^^ in the nucleus B, and Q^g in the nucleus A. As 
a result, the low- 2; gluon population is enriched in both nuclei. 



The next important observation is that the pr-broadening on such "ex- 
cited", or boosted nucleons, A^ is larger than in pA collisions, Apf,|^ > 
Ap'^lj^, since the density of target gluons is increased at small x. This should 
lead to a further mutual enhancement of broadening, i.e. a further increase 
of the saturation scales in both nuclei. Intuitively this seems to be clear, but 
a formal consideration below also supports this conclusion. 

Broadening is predominantly a process based on many soft rescatterings 
of the projectile parton. It was found in 15|, |l^ (see also [J]) that quark 
broadening is related to the dipole cross section, 

d(T{r, E) 



Qi^ = Ap'^E) = 2- 



dr"^ 



dzpA{b,z), 



(2) 



where pA{b,z) is the nuclear density at impact parameter b and longitudinal 
coordinate z, Since the process is soft, Bjorken x is not a proper variable, 
but instead the parton energy E should be used. The energy dependent qq 
dipole cross section was parametrized in the saturated form and fitted to 
photoab sorp tion, low Q'^ DIS data, and ttp total cross section in |28[ (see 



also |15|, |2J]). With that parametrization [15 
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where (T^^(-E) is the np total cross section; {tI^^t^ ^ 0.44 fm^ is the mean pion 
charge radius squared; QgN^E) = 0.19 GeV x {E/lGeV)'^^^* is the proton 
saturation momentum. 

Notice, that Eq. (|3]) Cg{E) has no scale dependence. It corresponds to the 
dipole-nucleon cross section crdip{r) = C^r^ and a soft scale characterizing 
the proton is implicitly contained in Cg. Strictly speaking, however, this 
coefficient is divergent at r — )■ 0, since it contains In(rr) |27|. This divergency 
originates from the ultraviolet behavior of the unintegrated gluon density 
J-'(x, kx) oc 1/k^ at large kx- In reality this divergency is not harmful due to 



the natural cut-offs discussed in [15|, and to a low sensitivity to their values. 
One should fix the In r dependent factor term in Cg at some value of r typical 
for the process under consideration, as is done in Eq. ([3]). 

The function Cg{E) in (|3]) describes broadening in pA collisions, result- 
ing from multiple interactions with "normal" target nucleons, whose parton 
distributions are the same as in pp interactions. However, as was discussed 
above, the wave functions of participant nucleons in AA collisions are biased 
towards a larger scale Ql and a higher multiplicity of the constituent par- 
tons. This is why the Cg in ([3]) acquires the second variable, a scale Q^ = Q^, 
Cg{E) =^ Cg{E,Q'^). To expose this scale dependence explicitly, we present 



the function Cq{E, Q"^) in the form [25|, l26 



vr 



Cg{E,Q^) = —a,{Q^)xgM{x,Q^) 



(4) 



where 5'Ar(x, Q"^) is the gluon distribution function in the target nucleon, and 
X = Q^/2mNE. 

One should be careful bridging Eq. (jl]) with its soft limit Eq. ([3]), since 
tts(Q^) and gNi.x^Q'^) are ill defined at small Q^. To regularize this problem 
we replace Q^ ^ Q^ + Qq, where Ql should be adjusted to the reproduction 



of the correct infra-red limit, 

— ttsiQl) xgN{x, Ql) = Cg{E = Ql/2mNx). (5) 

Apparently, the value of Q^ is not universal. It depends on s, x, and most of 
all on the PDF-analysis dependent behavior of gN{x, Ql), especially at small 
Q^. Here we follow this procedure. 

Thus, a participating nucleon simultaneously plays the roles of a beam 
and of a target. As a beam hadron its PDF is boosted to a higher scale 
due to multiple interactions it undergoes in another nuclei. As a target 
such a nucleon, being boosted to a higher scale Ql ^ Ql + Q%, it increases 
broadening of partons from another nucleus, since the factor Cg{E, Q^) Eq.(|l]) 
rises. This leads to a mutual enhancement of the saturation scales in both 
nuclei. Indeed, multiple rescatterings of nucleons from the nucleus A on the 
boosted nucleons in B proceed with a larger cross section, so broadening, i.e. 
the saturation scale in B increases, Q^^ ^ Ql^ > Q^^- For this reason, the 
nucleon PDFs in A get boosted more. Then the partons from B experience 
even stronger multiple interactions with such double-boosted nucleons in A. 
This results in an additional boost of the saturation scale in A, then, as a 
result, in B, and so on. 

Such a multi-iteration mutual boosting of the saturation scales is illus- 
trated pictorially in Fig. [3l where two raws of nucleons, Ta and Tg are dis- 
played on horizontal and vertical axes. In the left picture a gluon originated 
from one of the nucleons in A (the rightmost vertical dashed line) propagates 
through B interacting multiply via gluon exchanges (horizontal dashed lines) 
and increasing its transverse momentum squared by Ap"^ = Q^^. In a simi- 
lar way a gluon from B (the upper red line) interacts multiply propagating 
through A and gets broadening Q^^- ^11 multiple interactions occur at the 
soft scale Q^. The two gluons collide (the cross in the upper right corner) at 
hard scale Q^, but the original gluon distributions in the colliding nucleons 
are shifted to higher scales, Q^ + Q^^ ^^^ Q^ + Q% respectively. 

Then we observe that the t-channel gluon exchanges (vertical lines) be- 
tween the gluon originated from B and nucleons bound in A, become s- 
channel gluons propagating through B, after a Lorentz boost between the 
rest frames of nuclei A and B. Therefore, such gluons should also undergo 
multiple collisions with nucleons in B. This is shown by additional horisonal 
gluon lines in the right picture of Fig. |3l And vice versa, the original gluonic 
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Figure 3: Collision of two one-dimensional rows of nucleons Ta and Tb displayed on hori- 
zontal and vertical axes. Multiple interactions of colliding gluons propagated through both 
nuclei, including additional multiple scatterings of gluons, which carry out the interactions, 
are shown as is described in text. 



exchanges carrying out multiple interactions in the rest frame of B, propa- 
gate through A and also experience new multiple interactions, as is depicted 
in the right hand side of Fig. Ei These additional interactions boost each 
of the multiple interactions to a new scale, as well as the scale of the hard 
reaction Ql^^^^ => QIa{b)- 

In the next iteration (not shown in Fig. ^ the new exchanged gluons also 
experience multiple scatterings, and so on. The final gluon saturation scales 
Ql in the collision of two rows of nucleons Ta and Tb can be found solving 
the reciprocity equations {Cq =^ Cg = |Cg), 



QIb{xb) 



37r2 

2 

37r2 



"s(<5L + Qo) XBQNixB, QIa + Qo) Tb; 
^siQls + Qo) xa9n{xa, QIb + Ql) Ta- 



(6) 



These equations are the main result of the paper. Compared to Eq. ([3]) 
these equations take into account the modification of the properties of bound 
nucleons in each of the colliding nuclei due to multiple interactions in another 
nucleus and the following increase of the scale. Notice that the values of Qo 
in the first and the second equations ([6]) depend on xa and xb respectively, 
and therefore might be slightly different, but we use the same symbols to 
simplify the notations. 



The reciprocity equations should be solved numerically, but here we esti- 
mate the magnitude of the effect for the case of central collision of identical 
nuclei, i.e. Ta = Tb, and a glue-glue collision at mid-rapidity, xa = xb = x. 
In this case the system of equations ([6]) reduces to a single one, 

QL(^) = ^- «.(QL + Ql) x9n{x, QIa + Qo) ^A. (7) 

The scale characterizing multiple interactions of gluons is the mean trans- 



verse momentum of gluons (^t) ~ 0.65 GeV [28|, |29|. Therefore, the gluon 
distribution should be taken at x = {kx) / ^/s, which gives x = 3.25 x 10~^ 
and X = 1.18 x 10^^ corresponding to v^ = 200 GeV (RHIC) and 5.5 TeV 
(LHC) respectively. 

To proceed further we should fix the infra-red cutoff Qq given by Eq. (|5]). 
At the energy of RHIC the parameter in the right-hand side of Eq. (j5]) C{E = 
65 GeV) = 3.2. Then Eq. © results in Qli^/s = 200 GeV) = 1.84 GeV^. 
At the energy of LHC the factor C{E = 1.787 TeV) = 7.03. So we found 
Ql{y/s = 5.5 TeV) = 1.7 GeV^. These figures confirm our expectation of a 
weak energy dependence of the infra-red cutoff Qo- 

With these values of x and Qo we solved the equation ([7]) for a central 
collision of identical nuclei, relative to the modified value of the saturation 
momentum Qsa as function of Ta = Tb using the LO gluon distributions of 
the recent analysis MSTW2008 [23[ . The results are plotted in the left upper 
panel of Fig. H] as function of nuclear thickness Ta at the energies of RHIC 
and LHC. We see that the saturation scale of heavy nuclei may be as large 
as about 10 GeV^ at the LHC. 

To see the magnitude of the boosting effect we also show in the left bottom 
panel the boosting factor as function of Ta = Tb at the energies of RHIC 
and LHC. The enhancement is significant, especially at the energy of LHC, 
where it reaches a factor of three. 

We also solved the reciprocity equations ([6]) for non-central collisions, i.e. 
for Ta 7^ Tb- The boosting factor is plotted in the right panel of Fig. [Has 
function of Tb for fixed values of Ta = 2, 1, 0.5 fm~ from bottom to upper 
curves respectively. Solid and dashed curves show the boosting factor for nu- 
clei B and A respectively. These results confirm the expectation illustrated 
pictorially in the left panel of Fig. [2l Namely, in the limit Tg — )■ 0, corre- 
sponding to pA collisions, the boosting effect in nucleus A vanishes, while 
the parton distribution in the projectile proton is drifting to a higher scale, 
i.e. the boosting factor exceeds unity. At larger values of Tg the numerical 
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Figure 4: Left upper panel: the boosted values of the saturation momentum squared 
Ql^ calculated for T4 = Tb with Eq. ([7]) at the energies of RHIC and LHC as function 
of nuclear thickness. Left bottom panel: the boosting factors Ql^/Q^j^ as function of 
nuclear thickness. Right panel: The boosting factors as function of Tg for fixed values of 
Ta — 2, 1, 0.5 fm~ from bottom to upper curves respectively. Solid and dashed curves 
correspond to the boosting factors for nuclei B and A respectively. 

results in Fig. H] show that the saturation scales in both nuclei are boosted 
to higher values in accordance with Fig. [2l right. 



4. Gluon shadowing 

The gluon density at small x in nuclei is expected to be reduced com- 
pared with free nucleons. This phenomenon called gluon shadowing, is a 
part of saturation, but it also affects and diminishes the saturation scale. 
With this observation we formulated an equation for the modified saturation 
scale in pA collisions [4] . The found reduction of the saturation momentum 
is significant and is similar to the result of numerical solutions [22| of the 
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Thus, both Q'l^ and Q^^ plotted in 
m would be reduced in accordance with Ref. [^ (see 



Balitsky-Kovchegov equation 
the upper panel of Fig. 
Fig. 3 of that paper). 

However, the modification of the saturation scale due to mutual boosting 
in heavy ion collisions turns out to be practically unaffected by gluon shad- 
owing. Indeed, either the usual saturation scale Eq. (jl]), or the reciprocity 



equations ([SD are controlled by the gluon distribution function, which should 
be modified in a nucleus by a factor i?g(x, Q^,Ta), which is the nucleus-to- 
nucleon ratio of the gluon PDFs. The difference between the two is in the 
scales for gluons shadowing, which are Q^^ and Q^^ respectively. However, 
the scale dependence of as{Q^) g{x,Q'^) at small x, given by the DGLAP 
evolution, is very slow and can be neglected. Thus, the boosting factors 
depicted in the bottom panel of Fig. H] are not affected by gluon shadowing. 

5. Experimental observables 

Increase of the saturation scales in colliding nuclei should certainly lead 
to a rise of transverse momenta of produced hadrons in comparison with an 
extrapolation from pA collisions. However, final state interactions with the 
dense medium created in heavy ion collisions, significantly modify the trans- 
verse momentum distribution of produced hadrons. Exclusion is production 
of heavy quarkonia. Propagating through a dense medium they experience 
no energy loss, and their survival probability is practically independent of 
Pt- Therefore any observed modifications of the p-r-distribution of heavy 
quarkonia produced in nuclear collisions should be associated with initial 
state interactions. This fact makes them an excellent for study of gluon 
broadening in colliding nuclei. 

We expect an increased magnitude of broadening J/\l/ and T produced in 
AA compared with pA collisions for the same path length in nuclear medium. 
This could be interpreted as a signal of the effect of boosted saturation scale 
discussed here. 



Fig. E presents RHIC data |30| at ^i = 200 GeV for the mean J/^ 
transverse momentum squared versus the mean nuclear thickness covered by 
the projectile gluon in the rest frame of each of the colliding nuclei, calculated 
at impact parameter of collisions, b, corresponding to the measured centrality 
bin, 

{Ta + Tb) = -^ I dh [Ta(s)+Tb(6-s)1 Ta(s)Tb(6-s), (8) 



TAB{h) 

where TAsip) = J (fsTA^s) TBip—s). Data are compared with the prediction 
based on Eq. ([2]) applied to J /'^ production (see details in j^) depicted by 
solid and dotted curves for Au-Au and d-Au collisions respectively. These 
curves must coincide at ?/ = (left panel), but are different at forward 
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Figure 5: Data [30| for the mean transverse momentum squared of J/^ produced at 
y/s = 200 GeV an at y = (left panel) and y — 1.2 — 2.2 (middle panel) plotted as 
function of {Ta + Tb) defined in (|S]). Solid and dashed curves are calculated for AA 
collision without and with the boosting effect. Same for y^ = 5.5 TcV is shown in the 
right panel. The d-Au point in the middle plot should be compared with the dotted curve. 



rapidities (middle panel), where xai^xb, and we fixed rapidity at {y) = 1.7. 
We predict here only broadening, i.e. the slope of the curves, while the 
absolute value of {p^) , which is model dependent and has not been calculated 
so far, is treated as a fitting parameter. 

We also solved Eqs. [6] applied to J/\l/ production, and plotted the result 
by dashed curves. While at the mid rapidity (left panel) the effect boosting is 
sizable, although not strong, at forward rapidities (middle panel) it is hardly 
visible. This is a result of compensation between the effects of rising and 
decreasing Bjorken x is the two nuclei. 

Within the rather large error bars data agree with the theoretical expec- 
tations, but cannot resolve the weak boosting effect. Extrapolated to small 
values of {Ta+Tb) all curves should meet at the value of {p\) for pp collisions. 
We found {p'^)pp — 3-52 ± 0.3 and 3.06 ± 0.3 at |/ = and 1.7 respectively. 



which agree within the errors with the measured values [30|. 

A much stronger boosting effect for broadening of J /'^ is expected at 
LHC. An example at i/i = 5.5 TeV at ?/ = is plotted in the right panel of 
Fig. |5l Broadening in AA collisions is enhanced up to factor three compared 
with pA collisions at the same path length in nuclear matter. 

Notice, that this kind of enhancement for broadening in AA compared 
to pA collisions was observed recently in high statistics measurement of J/^ 
production in the NA60 and NA50 experiments at Eiab = 158 GeV |3l|- The 
magnitude of broadening in nuclear collisions was found twice as big as in pA 
measurements for the same path length in a nuclear medium. The magnitude 
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of the observed boosting is much larger than follows from equations (jS]) at this 
energy, and is probably related to another mechanism enhancing broadening 
in AA collisions due to interaction with gluons radiated in the preceding 



multiple collisions [32| . This mechanism correctly predicted the magnitude of 
the effect observed in the NA50/60 experiments. However, this contribution 
steeply falls with energy |32| , and is negligible at the energies of RHIC and 
LHC. 

Another observable sensitive to the saturation scale is hadron multiplicity 
J33| . In this case the boosting effect should lead to a jump of multiplicity 
in pA and AA collisions at the same number of participants. Indeed, such a 



discontinuity was observed [3J, |35[ in data for Au-Au and Cu-Cu collisions 
at a/s = 200 GeV in comparison with multiplicity dA and pp collisions. The 
magnitude of the observed enhancement is in accord with the boosting factors 
presented in Fig. HJ A detailed analysis of the data and comparison with 
theoretical expectations goes beyond the scope of this Letter, and will be 
published separately. 

6. Summary 

Nuclear targets have a larger resolution than a proton for parton distribu- 
tion in the projectile hadrons. As a result, the projectile parton distribution 
distribution is suppressed at large x and enhanced at small x. 

The increase of projectile parton densities in the case of nuclear collisions 
becomes a source of enhancement of the broadening experienced by the target 
partons propagating through the projectile nucleus, i.e. to an increase of the 
saturation momentum in the beam (see Fig. H]). Such a mutual enhancement 
of the saturation scales leads to the system of reciprocal equations (|6]). 

We solved the equations for central heavy ion collision and found the 
saturation scale for gluon radiation to be boosted up to a factor 1.5 at the 
energy of RHIC and a factor 3 at LHC. 

The saturation scale in nuclear collisions can be probed by measuring 
broadening of heavy quarkonia, which is not affected by final state interaction 
with the created dense medium. At the energy of RHIC the boosting effect is 
too weak to be observed in the currently available data for J/^ production, 
however it should be easily detected at the energies of LHC. 
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